Assessing reproducibility of inherited variants detected with short-read whole genome sequencing.

BACKGROUND: Reproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS. RESULTS: To dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×. CONCLUSIONS: Our findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

Diagnosis of genetic diseases in seriously ill children by rapid whole-genome sequencing and automated phenotyping and interpretation.

By informing timely targeted treatments, rapid whole-genome sequencing can improve the outcomes of seriously ill children with genetic diseases, particularly infants in neonatal and pediatric intensive care units (ICUs). The need for highly qualified professionals to decipher results, however, precludes widespread implementation. We describe a platform for population-scale, provisional diagnosis of genetic diseases with automated phenotyping and interpretation. Genome sequencing was expedited by bead-based genome library preparation directly from blood samples and sequencing of paired 100-nt reads in 15.5 hours. Clinical natural language processing (CNLP) automatically extracted children's deep phenomes from electronic health records with 80% precision and 93% recall. In 101 children with 105 genetic diseases, a mean of 4.3 CNLP-extracted phenotypic features matched the expected phenotypic features of those diseases, compared with a match of 0.9 phenotypic features used in manual interpretation. We automated provisional diagnosis by combining the ranking of the similarity of a patient's CNLP phenome with respect to the expected phenotypic features of all genetic diseases, together with the ranking of the pathogenicity of all of the patient's genomic variants. Automated, retrospective diagnoses concurred well with expert manual interpretation (97% recall and 99% precision in 95 children with 97 genetic diseases). Prospectively, our platform correctly diagnosed three of seven seriously ill ICU infants (100% precision and recall) with a mean time saving of 22:19 hours. In each case, the diagnosis affected treatment. Genome sequencing with automated phenotyping and interpretation in a median of 20:10 hours may increase adoption in ICUs and, thereby, timely implementation of precise treatments.

Bead-linked transposomes enable a normalization-free workflow for NGS library preparation.

BACKGROUND: Transposome-based technologies have enabled the streamlined production of sequencer-ready DNA libraries; however, current methods are highly sensitive to the amount and quality of input nucleic acid. RESULTS: We describe a new library preparation technology (Nextera DNA Flex) that utilizes a known concentration of transposomes conjugated directly to beads to bind a fixed amount of DNA, and enables direct input of blood and saliva using an integrated extraction protocol. We further report results from libraries generated outside the standard parameters of the workflow, highlighting novel applications for Nextera DNA Flex, including human genome builds and variant calling from below 1 ng DNA input, customization of insert size, and preparation of libraries from short fragments and severely degraded FFPE samples. Using this bead-linked library preparation method, library yield saturation was observed at an input amount of 100 ng. Preparation of libraries from a range of species with varying GC levels demonstrated uniform coverage of small genomes. For large and complex genomes, coverage across the genome, including difficult regions, was improved compared with other library preparation methods. Libraries were successfully generated from amplicons of varying sizes (from 50 bp to 11 kb), however, a decrease in efficiency was observed for amplicons smaller than 250 bp. This library preparation method was also compatible with poor-quality DNA samples, with sequenceable libraries prepared from formalin-fixed paraffin-embedded samples with varying levels of degradation. CONCLUSIONS: In contrast to solution-based library preparation, this bead-based technology produces a normalized, sequencing-ready library for a wide range of DNA input types and amounts, largely obviating the need for DNA quantitation. The robustness of this bead-based library preparation kit and flexibility of input DNA facilitates application across a wide range of fields.

Connexin Hemichannels: Methods for Dye Uptake and Leakage.

It is now clear that connexin-based, gap junction "hemichannels" in an undocked state are capable of opening and connecting cytoplasm to the extracellular milieu. Varied studies also suggest that such channel activity plays a vital role in diverse cell processes and abnormal hemichannel activity contributes to pathogenesis. To pursue fundamental questions in this area, investigators require methods for studying hemichannel permeability and dynamics that are quantitative, sensitive, versatile, and available to most cellular and molecular laboratories. Here we first provide a theoretical background for this work, including the role of cellular membrane potentials. We then describe in detail our computer-assisted methods for both dye uptake and leakage along with illustrative results from different cell systems. A key feature of our protocol is the inclusion of a mechanical stimulation step. We describe dye uptake, interpreted as connexin dependent, that is shown to be enhanced with reduced extracellular Ca2+, mechanically responsive, inhibited by TPA, inhibited by EL186 antibodies for Cx43 and sustained for more than 15 min following mechanical stimulation. We describe dye leakage that displays these same properties, with estimates of hemichannel numbers per cell being derived from leakage rates. We also describe dye uptake that is shown to be unaffected by a reduction in external Ca2+, insensitive to EL186 antibodies and relatively short-lived following mechanical stimulation; this uptake may occur via pannexin 1 channels expressed in the cells studied here. It is unlikely that cell damage plays a significant role in dye uptake following mechanical stimulation, given compelling results from various control experiments.

Whole-genome characterization of human and simian immunodeficiency virus intrahost diversity by ultradeep pyrosequencing.

Rapid evolution and high intrahost sequence diversity are hallmarks of human and simian immunodeficiency virus (HIV/SIV) infection. Minor viral variants have important implications for drug resistance, receptor tropism, and immune evasion. Here, we used ultradeep pyrosequencing to sequence complete HIV/SIV genomes, detecting variants present at a frequency as low as 1%. This approach provides a more complete characterization of the viral population than is possible with conventional methods, revealing low-level drug resistance and detecting previously hidden changes in the viral population. While this work applies pyrosequencing to immunodeficiency viruses, this approach could be applied to virtually any viral pathogen.

Gap junctions assemble in the presence of cytoskeletal inhibitors, but enhanced assembly requires microtubules.

The role of cytoskeletal elements in gap junction (GJ) assembly has been studied using Novikoff hepatoma cells treated with cytochalasin B (CB) to disrupt actin filaments or with colchicine or nocodazole to disrupt microtubules. After 60 min of cell reaggregation, freeze-fracture was used to evaluate quantitatively the "initiation," "maturation," and "growth" phases of GJ assembly. The development of junctional permeability to fluorescent dyes was also analyzed. The only effects of CB on the structure or permeability of the developing junctions involved an elongation of GJ aggregates and a small decrease in formation plaque areas. Colchicine (but not the inactive form, lumicolchicine) prevented the enhancement of GJ growth by cholesterol, but its effect on basal growth was equivocal. Nocodazole inhibited the growth of GJ, even under basal conditions, without an effect on initiation. Nocodazole also blocked the forskolin-enhanced increase in the growth of GJs and, in living MDCK cells, reduced the movement of transport intermediates containing green fluorescent protein-tagged connexin43. Thus, neither actin filaments nor microtubules appear to restrict GJ assembly by anchoring intramembrane GJ proteins, nor are they absolutely required for functional GJs to form. However, microtubules are necessary for enhanced GJ growth and likely for facilitating connexin trafficking under basal conditions.